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ABSTRACT 

In a case study (June 6-7, 2008) we report on how the internal structure of a coronal mass 
ejection (CME) at 1 AU can be anticipated from remote observations of white-light images of the 
heliosphere. Favorable circumstances are the absence of fast equatorial solar wind streams and 
a low CME velocity which allow us to relate the imaging and in-situ data in a straightforward 
way. The STEREO-B spacecraft encountered typical signatures of a magnetic flux rope inside 
an interplanetary CME (ICME) whose axis was inclined at 45° to the solar equatorial plane. 
Various CME direction-finding techniques yield consistent results to within 15° . Further, remote 
images from STEREO-A show that (1) the CME is unambiguously connected to the ICME and 
can be tracked all the way to 1 AU, (2) the particular arc-like morphology of the CME points 
to an inclined axis, and (3) the three-part structure of the CME may be plausibly related to 
the in situ data. This is a first step in predicting both the direction of travel and the internal 
structure of CMEs from complete remote observations between the Sun and 1 AU, which is one 
of the main requirements for forecasting the geo-effectiveness of CMEs. 



Subject headings: 
medium 

Introduction 



Sun: coronal mass ejections (CMEs) — solar-terrestrial relations — interplanetary 



Coronal mass ejections (CMEs) are violent ex- 
pulsions of plasma and magnetic flux from the 
outer solar atmosphere. Understanding their 
propagation characteristics and their internal 
structure is one of the main goals of the NASA 
twin-spacecraft Solar Terrestrial Relations Obser- 
vatory (STEREO, [Kaiser et aLl|2008[ ). Recently, 
there have been many efforts to obtain the direc- 
tion of propagation of CMEs both close to the 



2008} [Thernisien et al.|[2009l [Temmer et ar][2009] 
Davies et al.||2009 ). Using interplanetary observa- 



tions by the STEREO Heliospheric Imagers (His, 
Eyles et al.|2009| [Harrison et al.|2009| , [Davis et al. 



(2009) were able to predict the CME direction of 



Sun and in interplanetary space (e.g. Mierla et al 



propagation and arrival time at Earth with good 
accuracy. However, a major issue in forecasting 
the geo-effects of a given CME is its orientation 
and its internal magnetic field since long-duration 
southward fields lead to strong geomagnetic ac- 
tivity (e.g. 



Farrugia et al. 1993 Zhao & Hoek 



sema|1998| . Thus the above-mentioned techniques 
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need to be compared with the magnetic field and 
plasma signatures of CMEs as measured in situ 
in the solar wind (interplanetary coronal mass 
ejections, ICMEs). Those ICMEs containing a ro- 
tating magnetic field vector are particularly suited 
for these comparisons. These are generally called 
magnetic flux ropes (MFRs), and may be mag- 
netic clouds (MCs) if they satisfy other criteria, 
namely, above-average magnetic field strengths B, 
low proton temperature (Tp) and beta (/3p), and a 
large rotation of B (Burla ga et al.|[l"981 1. 

Using SOHO/LASCO coronagraph images, 
Cremades & Bothmer (20041 have shown that 



CMEs have a different morphology seen along or 
perpendicular to the axis of symmetry. In this 
Letter we relate, for the first time, both the mor- 
phology and direction of a CME seen with the 
STEREO-HI instrument imaging the heliosphere 
to properties of the associated ICME measured 
in-situ. Our study thus extends previous efforts 
to predict the magnetic field orientation in MFRs 
from solar disk and coronagraph observations (e.g. 



Bothmer fc Schw enn|1998||Yurchyshyn et al.|200T| 

Yurchyshyn, 2008) ) . The big advantage is that the 
CME can be seen in the HI images from a dis- 
tant spacecraft (STEREO-Ahead) passing over 
the other spacecraft (STEREO-Behind) at 1 AU, 
thus bridging the gap between remote and in-situ 
observations which existed before the STEREO 



era. 



The CME in question left the Sun - 21:00 UT, 
June 1, 2008 and its front boundary arrived at 
STEREO-B at 22:39 UT on June 6, 2008. In pre- 



vious work Robbrecht et al. (2009) used corona- 



graph and on-disk extreme ultraviolet (EUV) im- 
ages to show that this CME was not accompanied 
by typical on-disk signatures (filament eruption, 
flare, dimming, EUV wave), the main reason be- 
ing that the CME was of the streamer-blowout 
type lifting off slowly from high in the corona 
(i?0 « 1.15 — 1.4). This previously unrecognized 
type of CME is a good candidate for causing so- 
called "problem storms", i.e., those without an ob- 
vious solar origin (e.g. Zhang et al.||2007 ). In this 
Letter, we complement and extend the study of 



Robbrecht et all (2009 1 by modeling STEREO-B 



observations to obtain a more complete view of 
the MFR inside the ICME at 1 AU and discuss 
how its structure and direction could have been 
forecast from remote observations. CME-ICME 



events suitable for this kind of study are very rare 
during this unusually quiet solar minimum (see 
Davis et al.]|2009[ |Mostl et al.||2009a| ). This event 
is extraordinarily well suited for this kind of anal- 
ysis because (i) the absence of equatorial coronal 
holes (and thus fast solar wind streams) allows 
a straightforward interpretation of the HI images 
( Lugaz et al.|2008 1, and (ii) it is a slow CME, jus- 
tifying partly an assumption of constant velocity 
inherent to one of the applied techniques. 

2. In situ observations of the ICME at 
STEREO-B 

At 00:00 UT, June 7 2008, STEREO-B was 
25.38° east of Earth at 1.0545 AU and STEREO- 
A 29.26° west at 0.9567 AU, with a separation 
of 54.64°. Fig. [1] plots STEREO-B magnetic fleld 
data in Radial- Tangential-Normal (RTN) coordi- 
nates and plasma bulk parameters at 1 min res- 



olution, from the IMPACT/MAG (Acufia et al 



(2008); Luhmann et al. (2008 1 ) and PLASTIC 



(Calvin et al. 2008[ ) instruments, respectively. 



STEREO-B encountered a MFR between 22:39 
UT, June 6 and 12:27 UT, June 7, 2008 (inner 
two vertical solid lines) . This interval is partly an 
outcome of our reconstruction technique (see Hu 
et al.|[2004 l and partly determined by eye so that 
it encompasses the smoothly rotating B-vector. 

The two outer solid lines indicate a forward 
(ts=June 6 15:35 UT) and a very weak reverse 
shock tr= June 7 20:48 UT, see e.g. [Gosling et al.| 
( 1998 ) . There are two clear proton enhancements 
of Np « 30 cm-3 on either side of the MFR, the 
first one being, in part, the sheath behind the for- 
ward shock. The front of the second density peak 
arrives at tc ~ 12:00 UT, June 7. For most of the 
MFR, Tp, is higher than that expected for normal 



solar wind expansion (red trace, Lopez (1987)) 



This precludes the MFR from being a magnetic 
cloud. 

We modeled these data using Grad-Shafranov 



reconstruction (Hu & Sonnerup 20021 and force 



free fitting (the latter shown in parentheses, Lep- 
ping et al.|1990 1. (For recent multi-spacecraft val 



idation of the GS method see Liu et al.| ([2008|; 
Mostl et al.| ( |2009a|b| ).) The MFR is right-handed 
with an axis orientation of 6* = 51(37)°, (p = 
278(326)° in RTN coordinates (0 is the inchna- 
tion to the RT plane, (p is measured from R (0°) 
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towards T (90°)). Thus the MFR axis points 
roughly north-east as shown in Fig. [2] The mod- 
eled axial field strength Bo=l5A (20.3) nT, the 
radial scale size (diameter) D = 0.130(0.155) AU, 
and the impact parameter (the closest distance of 
the spacecraft to the MFR axis) isp = 0.81(0.60) x 
D. The toroidal and poloidal magnetic fluxes are 
$t = 0.72(0.37) X 10^1 Mx and $p = 1.19(1.47) x 
10^"'^ Mx/AU. These values are rather typical for 
ICMEs which are associated with medium-sized 
flares ( |Qiu et al.||2007 l. 

Figure [2] shows a 3D plot of the local MFR re- 
construction in the heliosphere, in a coordinate 
system centered on STEREO-B, seen from north 
looking down on the RT plane (close to the so- 
lar equatorial plane, top panel) and from another 
view point (bottom panel). Several fleld lines 
and its cross-section (color contour) are shown 
as well as its intersection (black contours) with 
the RT plane. It emerges from the reconstruc- 
tion that STEREO-B crossed the western flank 
of the MFR which intersected the RT plane be- 
tween 24° to 39° east of Earth, while the MFR axis 
crosses at 33°. This is important when compar- 
ing the ICME structure at 1 AU to the results of 
the CME direction-finding techniques in Section 
4. Further, in Fig. |2] the spacecraft are indicated 
as STEREO-B (blue triangle). Wind (black dia- 
mond), STEREO-A (red X). The Thomson sur- 
face for STEREO-A is the blue sphere, the yellow 
sphere the size of the C0R2 field of view (15 Rq). 
The CME directions from the method of triangula- 
tion (TR, yellow), forward modeling (FM, green), 
kinematic Fixed-<i> (KP, black) and elongation- 
fitting (EF, red) are indicated (see later sections 
and Table 1). The dashed (solid) black lines cen- 
tered on STEREO-A indicate the boundaries of 
the HIl (HI2) field of views. 



tion process, which includes corrections for flat- 
fielding, dark currents, background subtraction 



and the elimination of stars (Davies et al. 2009 



Brown et al. 2009 1 . Only then can intriguing de- 



tails of the morphology and evolution of the CME 
be seen. We strongly encourage the reader to look 
at the corrected running difference movies (see the 
electronic supplements), of which we show some 
selected images in Fig. [3] 

In HIIA and HI2A, both the CME leading edge 
(LE) and core clearly show an arc-like shape, typ- 
ical of a CME viewed orthogonal to its axis of 



symmetry (Cremades & Bothmer 20041, which is 



consistent with the derived MFR inclination of 
6* w 45° > 0°. The intensity of the CME ob- 
served in the HI white-light images is determined 



by Thomson scattering of electrons (e.g. Vourlidas 
[Oloward 2006) . Shortly after the CME enters 
the HI2A FoV it leaves the Thomson surface of 
strongest scattering (Fig. [2]) and the intensity de- 
clines. The northern part vanishes earlier (at ~50° 
elongation) than the southern part, which can be 
tracked all the way to STEREO-B. Assuming an 
evenly distributed density of the CME this effect 
is most likely caused by the CME and MFR incli- 
nation from which the angle tp might be inferred: 
the northern part is tilted away from the observer, 
looking to the east (west) this yields ip « 270(90)° 
for 6 « H-45°. This is again consistent with the 
MFR's moderate inchnation {9 « 45°). The as- 
sumption made here that STEREO-B crosses the 
CME close to its apex or nose is supported by 
direction-finding techniques in section 4. In this 
case (fi is expected to be close to 270 or 90 de- 
grees. Additionally, there is a 180° ambiguity in 
if depending on the sign of 9; for a southward- 
pointing axial field, the opposite is obtained (see 
also conclusions). 



Remote imaging of the CME by STEREO- 
A 



4. CME three—part structure and helio- 
spheric propagation 



Running difference images of the CME ob- 
tained from both cameras (HIIA: pointing 4 — 24° 
away from Sun center, 40 min cadence; HI2A, 
18.7 — 88.7°, 2 hours) of the Heliospheric Imagers 
aboard STEREO-A are shown in Fig. |3] from 
15:29 UT, June 2, when the CME completely en- 
ters the HIIA FoV until 16:09 UT, June 6 (w ts). 
The data have been put through the full reduc- 



We now compare the results of direction- 
finding results applied close to the solar surface 
(2 — 2Oi?0) with those used in the interplanetary 
medium (20 - 22Oi?0), and the ICME structure 
at 1 AU. 

First, the elongation e{t) of the LE and core 
were manually measured in the HI FoVs and fitted 
with the Sheeley et al. ( 2008[ ) formula (e.g. Rouil 
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lard et al.||2008| [Davies et al.|2009[ ) to infer a con- 
stant velocity and direction (see Table [l]) as well 
as the arrival times (18:34 UT, June 6 [LE] and 
13:08 UT, June 7 [core]). The technique assumes 
that the CME is well represented by the features 
tracked in the white-light images which result from 
line-of-sight integration and Thomson scattering 



effects ( Vourlidas & Howard 2006 1. The error bars 



for the direction and velocity follow from an error 
in e{t) of ±3°. The arrival times correspond very 
well with the arrival of the two density peaks on 
either side of the plasma void region (dashed ver- 
tical lines in Fig. [ij. It thus becomes plausible 
how, in this case, the classic CME features seen 
in coronagraphs ( jllUng fc Hundhausen||1985 1 cor- 
respond to the ICME observations. The extended 
bright front, with the LE as its outermost edge, 
corresponds to the interval of the first density en- 
hancement reaching up to the MFR; the void or 
dark cavity is the MFR; and the core is the dense 
material trailing the MFR. We speculate that each 
double density peak bracketing the MFR in Fig. [l] 
arises from material originating in the corona (in- 
nermost density peaks adjacent to the MFR) and 
solar wind swept up by the CME in interplanetary 
space (outermost peaks), i.e., the sheaths (see also 
Riley et al.|[2008| . 

Fig. |4] presents the time-distance plot from Sun 
to Earth. To convert e{t) to distance from the Sun 
r{t) in Fig. [4] two simple methods have been put 
forward: the "Point-P" and "Fixed-^>" methods 
dKahler & Webbl 120071 IWood et al.||2009||. Ba- 



sically, Point-P assumes that the CME's LE is a 
spherical front centered on the Sun, which yields 
( [Howard et allpOO? ! 



r(t) = d sin e(i) 



(1) 



where d is the distance of the observer from the 
Sun. Fixed-<i> assumes the CME to be point-like, 
propagating radially along a constant angle <!>' to 
the observer, 



rit) = -. 



sin (e + <i>') ' 



(2) 



The results are shown in Fig. |4] Linking the ob- 
servations with the timing tsi^c) when the CME 
LE (core) hits STEREO-B located at 226 i?©, the 
angle $ can be derived with $ = -30(-26)°. (An- 
gle 4> is measured from Earth while <!>' is measured 



from an observer, negative means east of Earth). 
In Fig. |2] the mean of the two directions is plot- 
ted. It is also seen in Fig. |4] that until r « 14Oi?0 
or 0.62 AU, the Point-P and Fixed-<I> methods 
yield quite similar values for r(i), with r(€) be- 
ing slightly lower for Point-P close to the Sun (see 
also Wood et al.||2009 |. After this, the unphysical 
assumption of Point-P that the CME cannot be 
tracked farther than d leads to a divergence from 
Fixed-(f> to smaller values of r(i). 

We now compare the directions obtained from 
the interplanetary techniques with direction- 
finding techniques derived from coronagraph ob- 
servations. Fig. |4] (inset) shows the kinematics of 
a distinct feature along the CME leading edge, 
one that could be observed and measured from 
LASCO C2/C3 and C0R2A images. As the fea- 
ture is seen from different viewing angles, the 
de-projected propagation direction of the CME 
is derived by applying a geometrical triangula- 



tion method as described in Temmer et al. (2009 1. 



The same feature has been tracked in the HIIA 
FoV (left yellow x symbols in Fig. |3]) which yields 
<I> = —45 ± 5°. The applied technique is lim- 
ited by the underlying assumptions (same CME 
feature distinguished in used instruments, radial 
outward motion of CME). However, the results 
are consistent with those derived by |Thernisien| 
et al. (20091, who use a forward modeling (FM) 



technique which presumes a particular shape of 
the CME, and find $ = -37± 10°. Table [l] sum- 
marizes the various results of CME directions and 
de-projected velocities, and the direction vectors 
are plotted in Fig. |2] In summary, the meth- 
ods yield roughly the same de-projected direction 
when compared to the MFR reconstruction, with 
a slight offset to the east by sa 10 — 15° of the 
methods applied closer to the Sun. Concerning 
the latitudinal direction, TR and FM yielded re- 
sults to within ±5° latitude of the solar equator, 
and the measurements for EF and KP were also 
taken within that range. This also means that 
STEREO-B encounters the apex or nose of the 
CME and not one of its "legs" . 

The CME velocities clearly show an acceler- 
ation between the coronagraphs and the hclio- 
spheric imagers as well as a consistency of the lat- 
ter with the proton bulk velocities as measured in 
situ. From EF a velocity difference of AV^ w 60 
km s~^ between the leading edge (401 km s^^) and 
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the core (340 km s~^) is consistent with AV = 49 
km s~^ between the MFR front and back bound- 
ary. These are values typical of ICMEs at 1 AU 
during solar minimum (AV ~ 60-70 km s~^: Jian 
et al.||2006|; (cf . 45 km s'^ found by [Burlaga fc 
Behannon 1982[ ). Note that both are larger than 
AV = 19 km s~^ between the averaged velocities 
of the two in situ density peaks which correspond 
to the white-light features. 

5. Conclusion 

We described observations of a CME as it prop- 
agates from the Sun to 1 AU and detailed the re- 
lationship to its modeled in-situ signatures. From 
this we found various interesting clues useful for 
forecasting the geoeffccts of CMEs: 

1 . How can the internal structure of the MFR 
be inferred from remote observations? The dis- 
tinct arc-like morphology points to an inclined 
MFR axis and, from the intensity evolution per- 
pendicular to the ecliptic plane (northern part 
vanishes earlier than the southern), we estimated 
« ±45°. The MFR's axial and poloidal field di- 
rections are consistent with the global solar dipole 



field (see the PFSS model in Robbrecht et al. 
2009) and indicated by black arrows in Fig. |2| 



The poloidal field goes from south to north and 
the eruption comes from the southern hemisphere, 
which yields a right-handed MFR (as observed), 
and thus an axial field to the east (away from the 
observer, thus 6 « +45°, ip « 270°). The MFR's 
axial field at the Sun is likely to be close to the 
solar equator. This can be seen in an accompany- 
ing movie in Robbrecht et al. (2009 1, which shows 



a more compact CME shape closer to the Sun in 
accordance with the 'Cremades & Bothmer ( 2004 1 



interpretation. A slight subsequent clockwise rota- 
tion (consistent with positive chirality for a twist- 
to-writhe conversion) of roughly 45° would then 
match the MFR observations. 

2. Is the direction of the CME obtained from 
various direction-finding techniques consistent 
with the observations at 1 AU? Yes, within the 
error bars all methods we use give a consistent 
picture. There seems to be a minor offset to the 
west when techniques are applied closer to the 
Sun. Whether this is due to a bias of the different 
techniques or a slight non-radial CME propaga- 
tion is the subject of further investigations. 



3. What is the magnetic fiux content of this 
pre-existing magnetic fiux rope on the Sun with- 
out reconncction signatures on the solar surface? 
We find that it is surprisingly high, comparable to 



events connected to C or even M-class fiares (Qiu 
et al.|2007[ ). Thus, this can be seen as evidence for 
a formation of the magnetic fiux rope before the 
eruption for this particular event. 

The tools and concepts we developed allow a 
first comparison between 3D CME directions and 
modeled in-situ data of those ICMEs which con- 
tain magnetic fiux ropes or magnetic clouds. Mod- 
eling ICMEs is important in order to obtain good 
estimates of their orientations and impact pa- 
rameters, and where the spacecraft intersects the 
axis (the eastern/western part of the central MFR 
south/north). Our results underline the need of 
having coronagraphs and instruments which im- 
age the heliosphere between the Sun and the Earth 
from a vantage view-point away from Earth so as 
to enhance our ability to forecast the geo-effects 
of coronal mass ejections. 
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Table 1: CME directions and velocities from vari- 
ous techniques. Angle <i> is the longitude measured 
from Earth (negative means eastward). Rq gives 
the distance from the Sun of the observations used 
for the respective technique and V the velocity for 
this distance range. 



Technique 


Instruments 


Rq 


<I>, deg 


V, km s-i 


Triangulation (TR) 


C2/C3 C0R2A 


3-14 


-45 ±5 


235 


Forward modeling (FM)'' 


COR2A/COR2B 


2-15 


-37 ±10 


265 


Elongation fitting (EF)'' 


HI1A/HI2A 


15-226 


-26 (-25) ±3 


401 (340) ± 15 


Kinematic Fixed-<i> (KP)^ 


HI1A/HI2A 


80-226 


-30 (-26) ±5 


440 ± 120^= 


In-situ modeling 


IMPACT/PLASTIC 


226 


-24 to -39 


403 / 392 / 384'! 



""taken from |Themisien et al.] ( |2009[ l. 

''Results for the CME core are in brackets after those for the leading edge. 

"^Here V is the median velocity from 80-226 Rq in a V{r) plot. In this range V(r) is roughly constant. 

''The in-situ V are means over the sheath region, the deHoffmann- Teller velocity of the MFR region, and a mean over the second 
density peak, respectively. 
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STEREO B 6/6/2008 22:39:0-6/7/2008 12:27:0 A t= 1 min 



20 




50 I 1 1 1 1 1 — I 1 n 1 1 1 r 




00:00 06:00 12:00 18:00 00:00 06:00 12:00 18:00 00:00 06:00 
Start Time 06-Jun-2008 00:00 UT 



Fig. 1. — Magnetic field and plasma data 
(STEREO-B). The interval between the two mid- 
dle solid lines is used for the reconstruction, the 
two outer solid lines indicate the forward and re- 
verse shock. Dashed lines are the arrival times 
from the elongation fitting method for the CME 
leading edge (left) and CME core (right). From 
top to bottom: magnetic field magnitude and 
magnetic field components in RTN coordinates (R 
pointing radially away from the Sun, T is the cross 
product of the solar rotation axis and R, N com- 
pletes the right-handed triad), proton bulk veloc- 
ity, proton number density, proton temperature 
(black) and expected temperature (red), proton 
beta and the total, magnetic and plasma pressure. 
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ICME 6-7 June 2008 




Fig. 2. — 3D plot of the local MFR reconstruction 
in the heliosphere (color contour) in a coordinate 
system centered on STEREO-B, seen from north 
looking down on the RT plane (close to the so- 
lar equatorial plane, top) and from another view 
point (bottom). For further details see text. This 
figure is also available as an mpeg animation in the 
electronic edition of the Astrophysical Journal. 
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HIIA 2008-06-02 15:29 | HI1A 2008-06-03 01:29 




Fig. 3.— Evolution of the CME in STEREO-A 
HIl (top 4 images) and HI2 (bottom 4 images), 
cf. Fig. |2] Earth (E), STEREO-B (B) are indi- 
cated as well as the elongation of Mercury. The 
tracked features of the CME leading edge and core 
for Fig. |4] are given by yellow crosses. The last im- 
age was taken approximately at the shock arrival 
time at STEREO-B (15:35 UT). Note the distinct 
appearance of the CME in both FoVs as an arc- 
like shape, i.e. a CME seen perpendicular to the 
axis of symmetry. This figure is also available as 
an mpeg animation in the electronic edition of the 
Astrophysical Journal. 
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Fig. 4. — Time-distance plot of the CME lead- 
ing edge and core. The elongation of the fea- 
tures in the HI images has been converted to dis- 
tance by using the Point-P and Fixed-<i> methods 
($ = -30° for the leading edge and $ = -26° for 
the core). The in-situ measured arrival times 
(shock) and tc (core) are indicated as 'x' at the 
position of STEREO-B (226 Rq). The inset shows 
the STERE0-A/C0R2 and SOHO/LASCO ob- 
servations used for the TR method. 
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